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ABSTRACT 



An automated digital computer technique of control 
system design is presented. The emphasis is on 
compensator design but the method is applicable to the 
design cf any system with free parameters. Signal 
representation and system response are in the time 
domain . 

Tte inputs required from the engineer are the 
system configuration, the desired output response and 
the free parameters. A parameter minimization routine 
is then used tc minimize a specific cost function and 
to set the free parameters. A graphical output of the 
desired response and actual system response is then 
produced fcr comparison by the engineer. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 10 

II. PROGRAM DEVELOPMENT AND IMPLEMENTATION 12 

A. GENERAL 12 

E. MAIN PROGRAM 15 

C. PLANT 15 

1. Block Data Card 16 

2- Block Connections 17 

3- Drives 17 

4- Standard Transfer Function Blocks 18 

5' Parameters to be Optimized 19 

D. DESIRED RESPONSE (XDATA) 20 

E. CCST FUNCTION (PERFORMANCE INDEX, PI) 20 

F. CINIMI ZAT ICN ROUTINE 22 

1. General 22 

2. Explicit/Implicit Constraints and 

Start Points 22 

G. GRAPHICAL OUTPUT 23 

III. INVESTIGATION OF PROGRAM PERFORMANCE 25 

A. DATA INPUT, AN EXAMPLE 25 

B. TACHOMETER FEEDBACK 28 

C. CASCADE COMPENSATION 34 

D. CASCADE LEAD COMPENSATION 41 

E. SERVO SYSTEM COMPENSATION 47 

IV. DISCUSSION, CONCLUSIONS AND RECOMMENDATIONS ... 61 

A. DISCUSSION 61 

B. CONCLUSIONS 64 

C. RECOMMENDATIONS FOR FUTURE WORK 65 

V. APPENDIX A BLOCK DATA CARD FORMAT 67 

VI. AEPFNDIX E PLANT FLOW CHART 69 

VII. APPENDIX C CADS MODIFICATIONS FOR 



5 



PBCEIEH III. E 80 

VIII. CflIS EROGEAH LISTING 82 



6 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 



14 

24 

26 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

33 

40 

41 

42 

43 

44 

45 

46 

48 

49 

51 

52 



LIST OF FIGORES 



Typical System to be Optimized 

CADS Erogram Flow Chart 

fiard-Leonard Speed Control System 

lachcmeter Compensated System 

CADS Elcck Diagram of Tachometer System 

CADS Compensated Tachometer Feedback System 

Response 

Type Three Block used for Pseudo Tachometer 

Feedback 

Pseudo Tachometer Feedback 

Eseudc Tachometer Compensated System's 

Response 

Third Order Plant to be Compensated 

lag Compensated System 

Uncompensated, Compensated BODE PLOT 

Conventionally Compensated System Response .... 

CADS Lag Compensated System 

CADS Compensated System Response 

Plant for lead Compensation 

Conventional Lead Compensated System Response 

CADS Elock Diagram for Lead Compensation 

CADS Lead Compensated System Response 

CADS Compensated System Block Diagram 

CAES Relaxed Boundary System Response 

Serve Drive Mechanism, Original Compensation 

CAES Eagram of Servo Drive Mechanism 

Original Servo System Response 

CADS Compensated C System Response 



7 



26 CAES Compensated C -1 System Response 54 

27 CAES Compensated C -2 System Response 55 

N 

28 CADS Optimized System Response 57 

29 Serve System Response, CADS Values with 

a = 0.5 60 

30 Elcck Reduced Ward-Leonard Drive 63 



8 



ACKNOWLEDGEMENTS 



The author wishes to thank Dr. George Thaler for his 
guidance and encouragement during the course of this 
investigation. Kris Butler, Ed Donnellan and M. Anderson of 
the W. B. Church Computer Center deserve special recogrition 
for their patient and professional assistance. Lastly, I 
thank my wife, Patricia, without whom this thesis could 
never have teen accomplished. 



9 



introduction 



I. 



The past two decades have been witness to an ever 
increasing use cf the digital computer. Engineering usage 
and design are prctably the most important justification for 
the large memory, high speed computers of tcday. Classes in 
computer application to engineering problems have fceccme 
part of tfce established curricula at most universities. 
Numerous papers have appeared on the application of the 
computer to engineering problems £ 1 ] which until recently 
were unsclvable or at best required long, tedious 
procedures which gave only approximate solutions. 

Control system engineering has relied increasingly on 
computer simulation of large scale systems to verify that 
design specifications have teen met and tc make 
modifications to a system even before producing a prototype. 
There are programs available to simulate virtually any 
electrical circuit, or control system, either in transfer 
function form or as a system of first order differential 
equations. ethers draw the Bode, Nichols or Nyquist plots 
cf open and closed loop systems. Some programs help design 
compensators which use an iterative method to achieve the 
desired frequency response [2] [3]* Researchers continue to 
better adapt the computer to engineering usage. 

Cant alapiedra £4] has used an iterative method tc find 
the optimum reduced order model for large order systems. 
MacNamara [5] went even further and used an iterative method 
tc find the optimum compensation for an aircraft autopilot. 
It would appear that these techniques could he extended and 
applied to the direct simulation and design of control 
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system circuits in the time domain. 



The irtert of this thesis was to develop a user oriented 
program which cculd simulate a wide variety of control 
systems and determine the values of the gains, poles and 
zeros necessary to produce a desired response. A convenient 
means of data card input was to be provided to specify (1) 
the control system which was to be optimized and (2) the 
desired response. A locally available function minimization 
routine (ECXILX) was to be used to optimize the simulated 
system's output. 

To simulate the system which is to be optimized, 
transfer functions which are commonly encountered were 
reduced tc first order linear differential equations. These 
equations were then programmed so. that the transfer function 
blocks cculd be connected in an arbitrary fashion by data 
card input. Several common nonlinear transfer blocks were 
also provided. The program will simulate the system with 
the known parameters and then allow all unknown or 
adjustable parameters to be fixed by the computer 
optimization routine to achieve the desired response. 



1 1 



II . 



PRCGHAM DEVE LOPMENT AND IMPLE MENT ATION 



A. GENEPAI 



Any program which is to be of maximum tenifit tc the 
user must have a simple means of data input and an output 
which is- easy to interpret and apply to the problem at hand. 
The input data should have a physical significance that does 
net lose its relevance through the programming of numerous 
eguatiens. The program should be a readily usable tocl and 
not a problem in itself. The intent of this thesis is to 
'present such a program. Computer Automated Design of Systems 
(CADS) , which is readily used for simulation with 
optimization of the output. Optimum in the sense that the 
output response cf a given system configuration is as near 
the desired response as possible. 

Most ccntrol system design starts with a preposed 
transfer function block schematic to achieve a desirec time 
domain response. CADS has the common transfer functions 
built into the pregram and the input data can come directly 
from the proposed system schematic. The transfer blocks 
which are available for system simulation are presented in 
Table I. These blocks should be adequate, either separately 
or in various cascade combinations, to represent most 
ccntrol systems. 
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ll£I 

CF SYHBOL EQUATION SOLVED 

EICCK 
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To use the program the engineer must know the system 
configuration in transfer block form and the desired second 
order output response. If other than a second order 
response is desired, this may be easily specified but 
requires seme knowledge of how to input the desired 
response. CAIS will take the transfer block system, connect 
it, compare the given system response tc the desired 
response and set any free parameters to achieve the closest 
match possible of the system output to the desired response. 

Figure 1 is representative of the type of systen the 
pregram car simulate and optimize. In figure 1, the 
parameters of the numbered transfer function blocks are 
known or fixed by equipment limitations. Transfer blocks X, 
I, and Z contain variable parameters which must be selected 
by the program to make the system output reproduce a desired 
output function as closely as possible. 




FIGQBE 1 Typical System to be Optimized 
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flAIK FBCGFAM 



E. 



The KAIK program was developed to control the selection 
of the sutfurcticns used in the optimization process ana to 
compute tie cesired response of the systems which were to be 
optimized. A second order step response is calculated from 
eguation (I) as the desired response and stored in the array 
XLATA. 



% 



G 

u (t) 

S2 + 2 8 HrjS + 'W 2 



(I) 



The MAIN is essentially a bookkeeping routine which controls 
the program execution. The subfuncticn operations it 
controls are defined in the following sections. 



C. ELANT 



The ccmmcn simulation routines LISA, DSL, ECAP, CSEP and 
INTEG were investigated in an attempt to adapt them tc the 
general problem specified above. These programs reguired 
either the input cf the system's state variable equations, 
Laplace transform equations or the node-transfer function 
pair for each system simulated. These programs are useful 
for simulation if all of a system's parameters have been 
specified. Eowever, each simulation is problem specific and 
these techniques are not readily integrated with other 
subfuncticn programs. A simulation subprogram which was 
capable cf accepting data card input tc simulate a wide 
variety cf ccntrcl systems and of working with other 
existing subprograms had to be developed. 
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Subfuncticn PLANT was developed to simulate the systems 
which were tc be optimized. It provides the versatility 
necessary to simulate numerous system configurations and is 
capable cf working with a minimization routine to optimize 
the variable parameters. The program reads the transfer 
function block connections from data cards. This data is 
then used to automatically set up the system equations in 
state variable format. These state variable equations are 
solved as first order, ordinary differential equations by 
the Bunge-Kctta-Gill forth order method. The programming 
was done in fortran IV and all calculations are in double 
precision. The capability of connecting the transfer 
function blocks in any configuration (feed forward, 
f eedback ,etc .) is possible with this program. Simulation 
flexibility is provided by having each block within the 
system capable of accepting an external forcing function. 

1 . fleck Data Card 

Several possible ways of inputting the data 
necessary tc simulate arbitrary system configurations were 
investigated. A means of defining a system configuration 
directly from a block diagram schematic was chosen as the 
most preferable method. Using this approach, each data card 
was designed to be directly related to a transfer furction 
of the system. This resulted in having direct access to the 
transfer function parameters for optimization. 

Tc simulate the system, a data card is prepared for 
each transfer block in the schematic diagram. The data card 
certains a field cf numbers which specify the block number, 
the type cf transfer function contained within the block, 
the input node and the output node to which the block is 
connected and the values of any parameters associated with 
the transfer function. The general format of the data 
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cards used tc input the system configuration is shown telcw. 



|1 J11 |20 1 40 |60 (Column Number) 

£LKCCC = E E V V G P Z 

Where : 

CC = Position number of the block 
D = lype of block (number) 

EE = Input node number 
VV = Output ncde number 

G # E, and Z are parameters cf the transfer function. 

2. Elcck Conn ection s 

Ihe number of transfer blocks in the system and the 
data cards associated with each of the blocks are read upon 
initial entry into the simulation routine. The program then 
connects the transfer function blocks in the proper order by 
comparing the input node number of a block to the cutput 
ncde number cf every other block. Whenever these two 
numbers are egual, the blocks are known to be connected and 
a flag is set egual to 1.0. The flags are then used to 
identify the input drives to each block. 



3. Erives 



The input guantity to each block is called the 
EBIVE. The program was designed to allow for multiple 
inputs to the system being simulated. This flexibility was 
achieved by making the input to each block egual to the sum 
of the outputs cf all blocks connected to the input node 
plus any external forcing function (DBVIN) feeding the 
transfer blcck. The input DRIVE to a block is determined as 
shewn in equation (II) . 
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EBIVE (i) 



THA (j) *FLAG (j,i) + ERVIN (i) 



(II) 



THA(j) is the output of block j. FLAG ( j , i) is 1.0 if block 
j is connected to block i and zero if it is not connected. 
ERVIN (i) is any external forcing function specified by the 
user which drives block i. DRVINs must be inserted in 
subroutine fLANT as Fortran IV statements. The standard 
program has ERVIN (1) = 1.0 specified as a unit step input to 
the system. Problem III-E in the section Investigation of 
Program Eerfcrmance demonstrates how multiple, time varying 
inputs are tc be inserted in the program. 



4 . 



Standard Transfer Function Blocks 



The transfer function blocks available for system 
simulation are shown in Table I. These blocks were selected 
because of their common usage in the modeling process. They 
were also fcund to be adequate either separately or in 
various cascade combinations to represent most control 
systems. 



The transfer function equations for each tjpe of 
block were written in state variable format and stored in an 
array named THA. The program reads a number from a block's 
data card which specifies the type of transfer function 
associated with the block. The system equations are then 
solved sequentially by selecting the type of equation 
assccitaed with block number one, solving fcr its output and 
then sequencing tc block number two, etc. The integral 
equations are solved by a modified RKL fourth order method 
in the suhfunction routines. Three integration routines 
were necessary tc store the intermediate results obtained 
for those equations involving a double integration. The 
four sutf u notions , RKLDE2, RKLBE3 , CCPLX and RKIEE4 are 
called by ILAtiT. EKLDE2 is used for the integration of a 
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type two 
three blcc 
subf unction 
integration 



block. BKLDE3 is used for integration of a 
k and to store intermediate guantities 
CCPLX . CCPLX calls RKLDE3 and RKLDEh 
cf a type four block. 



type 

for 

for 



No prevision has been made for the input of initial 
conditions cr the integrators. Therefore, the integrations 
must start at time equal zero. The user must specify the 
step size to be used for integration (DT) and the problem 
run time (If) . To conserve computer time ET should be made 
as large as possible. DT = TF/1000 is suggested as 
appropriate in most cases. Equally important is that TF be 
as small as possible.. Ose of the program has shewn that 
letting TF be greater than the transient response time of 
the system is rarely justified. For preliminary analysis TF 
was kept to only slightly longer than the time of the first 
undershoot for second order dominant systems. 3ecause of 
the complexity of subroutine PLANT a flew chart of PLANT has 
been included as Appendix B. 



5 • Param eter s to be Op timi zed 



The parameters of the system which are to be 
optimized by the minimization routine EOXPLX must be 

specified in PLAN 1 . The minimization routine returns the 
trial values cf the variable parameters to PLANT in an array 
labeled 'C'. Regular Fortran IV statements equate the 
variable parameters to the values in this array. If a 
pcle-zerc pair were to be optimized, the following 
statements would be inserted in PLANT: P (i) = C ( 1 ) 

and Z (i) = C (2) . The location for the preceding statements 
is clearly indicated in the program listing for PLANT. 
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D. 



DESI BED RESPONSE (XDATA) 



The criteria against which the system response will be 

compared will vary according to the application cf the 

system. Since most design work on control systems is done 

cn the basis of second order dominance, the program was 

written tc provide a second order step response with 

adjustable 5 , W and gain as the basic criteria against 

n 

which the simulated system will be compared. The desired 6, 

H and gain are read in as input data. The step response is 
n 

then computed in the main program by subfunction RKLEEC and 
stored in the array called XDATA. Any ether time domain 
response nay be specified by the user by removing the second 
order step response equations and replacing them with the 
equations of the desired response. An example is provided 
by problem III-E in Investigation of Program Performance. 
If the pregram is being used for simulation only and nc data 
curve is desired, setting the input variable LEAP = 1 will 
cause the program to bypass the computation cf the standard 
second craer data equation. 



E. C0S3 FOKCTTON (PERFORMANCE INDEX, PI) 



The achieved system response is compared with the 
desired response and the difference is the error. The 
program searches for the parameter settings which will 
minimize this error. The cost function may be specified by 
the user tc weight the system outputs as desired in 
subfuncticn FE. The default cost function cf the program is 
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1 1 € integral error sguared. 
of a weighted cost function 



TF 



J = J (Err) 2 dt. An 
is presented in Section 



example 

III. 
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E. 



MINIMIZilIGN ROUTINE 



1 . General 



Ihe free parameters are optimized to reduce the cost 
function by the complex method of M . J. Box. [6] Box's 
constrained optimization method has been programmed at the 
Naval Postgraduate School by R. R. Hilleary as a subroutine 
called ECXELX. This subroutine will find the minimum of an 
arbitrary function (cost function) subject to arbitrary 
explicit constraints and for implicit constraints. Explicit 
constraints are defined as upper and lower bounds on the 
free parameters. Implicit constraints may be arbitrary 

functions of the free parameters (e.g. P P < 178). 

1 2 

Tvc function subprograms are used to evaluate the 
objective function and implicit constraints, FE and KE 
respectively. The method EOXPLX uses to search for the 
values of the free parameters which minimize the cost 
function is explained in the computer program listing. 

2 . Expiic it/Implicit Cons tr a ints an d Sta rt Po in ts 

ECXELX searches a feasibility region (n-dimensicnal 
space, where n = number of free parameters) defined by upper 
and lower bcunds on the free parameters for a minimum cost 
function. Ihe smaller the region defined by these 
boundaries, the more rapidly the program will converge to 
the optimum parameter settings. Good engineering judgement 
will be necessary to keep the feasibility region as small as 
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possible. The boundaries of the search region are read from 
data cards by the main program as the upper bound (XO) and 
the lower bcccd (XL) for each free parameter. 

Implicit constraints may be any arbitrary function 
of the free parameter desired. If an implicit constraint 
such as the product of a pole-zero pair must be less than 
seme number is to be evaluated, it must be supplied by the 
user to the subfunction KE. No implicit constraints were 
used in this thesis. 

The starting values cf the free parameters (XS) are 
read in by the HAIN frem data cards. A good choice of 
starting values will dramatically reduce the time required 
for optimization. A preliminary root locus or Bode plot 
method of estimating the best values of the free variables 
should be accomplished whenever possible. 



G. GRAPHIC AI OOTEUT 



Two subroutines were written to provide for the 
graphical output of the desired response and the best system 
response achieved by the optimization process. The user may 
select, by data card input, either subroutine PPL! which 
provides a high speed printer plot or subroutine PIC which 
provides a calcomp graph. Every fifth integration point 
stored in the arrays XDATA and THAOUT is plotted. 
Subroutines PELT and PIC call the subroutines PLOTP and CRAW 
respectively. PLOTP and DRAW are standard plotting routines 
at the NP£ computer facility and are not a part of the 
simulation program. Figure 2 diagrams the information flow 
and data input to the program. 
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Optimization 

Simulation 



FIGCSE 2 CADS Erograa Flow Ciart 
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III. 



INVESTIGATION OF PROGRAM PEEFORHANCE 



The example problems presented below were used tc aid in 
the development cf the optimization program. The order of 
difficulty cf the problems progresses from a simple text 
bock single variable, single input system tc a multivariable 
operational servo drive system which has multiple inputs and 
discrete level feedback. An example of how a schematic 
diagram representation of a system is prepared for input to 
the program is presented prior tc considering the example 
problems . 



A. EAT1 INECT, AN EXAMPLE 



The Rard-Ieonard drive system [7] shown schematically in 
Figure 3 (a) has two variable parameters. The gain cf the 
amplifier ana the tachometer feedback are available to 
adjust the system's response. To simulate the system a 
block diagram representation of the system is drawn as shown 
in Figure 3 (c) using the transfer blocks from Table I. 
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X 

(C) 



fIGCBE 3 Hard - Leonard Speed Control System 
Using Feedback. (A) Schematic Eiagram 
(B) Elock Diagram (C) Blcck Diagram 
Using CADS Blocks. 

Where G(1) = K/r P(1) = 1/ T 

G (2) = Km/Ea 

G (3) = 1/J P (3) =0 

G<4) = -Km 

G (5) = -Kt 

DEVIN (1) = Er 
DEVIN (3) = -T t 

Ihe nodes cf the block diagram are then numbered 

sequentially 1,2, The blocks between the nodes are 
also numbered sequentially 1,2,...N as shewn in 

Figure 3 (c) . Data cards (N) are then prepared for each 

blcck which specify the block number, type cf blcck, the 
input nede, the output node, and the parameters contained 
within the blcck. In Figure 3 (c) , for example, block 1 is 
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a type two transfer block connected between nodes 1 and 2. 
lie data c arc input for this block would be 

BIKC 12=0102 K/ t 1/ t 



The program reads the data card input and connects the 
blocks by setting a FLAG = 1. whenever the input node 
number to a block is the Same as the output ncde number of 
any ether block. The input to a transfer block is then 
determined tc be the sum cf the outputs of all blocks 
connected tc the input node plus any external forcing 
function driving the input, node. The program has a unit 
step specified for DEVIN (1). If this is the only input to 
the system, nc action is necessary on the part of the user. 
If ether than a unit step input to the system is desired or 
if there are other external forcing functions such as 
EEVIN (3) , they must be specified and placed within the body 
cf the subfurction PLANT as Fortran IV statements. For the 
example shewn in Figure 3 (c) , a card with the equation 

DEVIN (3) = ffTJ 



would have tc be inserted preceding the drive equations. An 
example cf hew multiple, time varying drives are specified 
is given in section III. E. 



fihen optimizing a system, seme of the input quantities 

will be unknown or variables. These variables must also be 

* 

assigned within subroutine PLANT. To optimize the variables 

K and K cf Figure 3 the following two statements would be 
1 t 

inserted in PLANT: 







G (1) 


= C(1) 






G (5) 


= C (2) 


where C ( 


1) 


and C (2) 


are the variables which will be 


optimized 


by 


subroutine 


BOXPLX . 
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E 



IflCBCMf 1£E £E£DBACK 



The first optimization problem attempted was one which 
had an exact solution that can be found by algebraic 
methods. An instrument servo [8] with unity feedback and 
fcrward transfer function 

1000 

G (S) = 

S(S+10) (III) 

was to be compensated with tachometer feedback as shewn in 

Figure 4. The only specification for the system's 

performance cf this single variable, seccnd order system was 

the simple requirement that the closed lccp roots have a 

5 = 0.7. 




£ IGUHE 4 Tachometer Compensated System 



The system shewn in Fig 
order tc achieve the 
characteristic equation of 

S 2 + (10 + 103 

The required K = 0.0343 
IV. t 



re 4 was redrawn as Figure 5 in 
tachometer feedback. The 

the system shown in figure 5 is 

K ) S + 103 = 0 (IV) 

t 

may be calculated from eguaticn 
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FIGOBE 5 CADS Block Diagram of Tachometer System 

CADS was programmed to optimize the system with the 
variable, , specified to be between the limits 

0.C1 < K < 1.0. The desired response was specified to be 
t 

the standard second order step response for 6 = 0.7, 

Wn = 100C. CADS determined the optimum value for to be 

K = Q.03h3. Figure 6 shows the system's step response and 
the desired response are virtually superimposed. 
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THAOUT 




EIGOfi 6 CADS Compensated Tachometer Feedback 

System Response 
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The type ooe system, in the preceding example, allowed 
derivative feedback from the forward path without reguiring 
a block which was capable of differentiation. One may 
encounter a type zero system or a system which cannot be 
configured to provide derivative feedback from the forward 
path. The possibility of providing a pseudo derivative 
feedback using a type three block was investigated for these 
cases . 



A type three block with Z = 0 and G = P is shewn in 
Figure 7. If the pole is placed far out on the real axis, 
this block approximates derivative feedback. 
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FIGUBE 7 Type Three Block used for 
Pseudo Tachometer Feedback 

The effect of using this pseudo tachometer feedback on 
the complex roots of the closed loop system is negligible. 
It does add c real root at P = -964. 

The system of Figure 5 was redrawn as shown in Figure 8 
using the pseudo tachometer feedback. 
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FIGURE 8 Pseudo Tachometer Feedback 

The optiiiizaticn program calculated K = 0.0352 fcr the 

t 

above configuration. This gave a 5 = 0.715. 



The system response and the desired response are shewn 
in figure S. The two responses are again rearly 
sup erimpesed . This method of providing derivative feedback 
has the disadvantage of adding an integration step tc the 
problem solution with the concomitant increase in problem 
solution time. 



32 



THAOUT 




Time (sec) 

fIGURE 9 Eseudo Tachometer Compensated System's Response 
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c. 



CASCADE CCMEEMSATION 



Having proven the feasibility of the program optimizing 
a single variable system where an exact solution was 
available, tie next problem considered extending the problem 
scope to a third order plant with two variables. The 
unstable plant that was to be compensated is shewn in 
Figure 10. 
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FIGOBE 10 Third Order Plant to be Compensated 

The plant was tc be stabilized using a single section 

cascade compensator. The compensated plant was required tc 

have H <2, without reducing the error coefficient. [8] The 
PH 

compensated system is shown in Figure 11. To keep the error 
coefficient constant, the compensator used was a simple lag 
network ( r,< r % ) . 
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FIGOBE 11 Lag Compensated System 

A compensator which would meet the recuired 
specif icaticrs was calculated by conventional methods as a 
check oe the program's performance. The Eode plot c£ the 
uncompensated and conventionally compensated system is shewn 
in Figure 12. The values of r ( and r. for Figure 11 which 
would meet the design requirements were determined to be 
t, = 1C., t z = 100. The compensated system's response 
using these values for the compensator is shewn in Eigure 
12 . 
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The conpensa ted system was redrawn as shewn in 
Figure 14 using the standard program transfer blocks 
available for system simulation. 




FIGURE 14 CADS Lag Compensated System 

A staraard second order response of <5 = .3, 

W = 0.75 was chosen as the desired response for the 
n 

optimization program to match. The limits placed on the 



optimization program were .001 < P < .1 and .01 < Z < 1. 

C C 

Arbitrary values to begin optimization were specified as 

E = .01 anc Z = . 1 in the logarithmic centers of the 
Co Co 

search 2 cnes. CAES determined T t = 8.078 and f 2 = 67.47 



as the optimum parameter settings. Figure 15 shows the 
desired response and the program compensated system 
re spense . 



The response of the system compensated by CADS is more 
nearly the desired response than is the conventionally 
compensated system. One should remember that the specified 
response is for a true second order system whereas, the 
compensated system is forth order. Therefore, a perfect 
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match cf desired versus actual responses could 
obtained . 



net be 



39 



THAOUT 




F1G0EE 15 CADS Compensated System Besponse 
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CASCADE LEAD COMPENSATION 



The complexity of the next problem 
was extended to five free parameter 
figure 16 was to be Used to foil 
sine-wave input of 200 rad/sec. The o 
be almost exactly the same as the inpu 
output cculd not lag the input b 
sections cf cascade compensation were 



to be solved by CADS 
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The plant 


shewn 


in 
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amplitude 
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amplitude was 
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amplitude. 
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FIGUBE 16 Plant for Lead Compensation 



The specified requirements may be interpreted as an open 

lccp gain > 15 db and a phase angle 90° at W = 200 from 

n 

a Nichols plot. A cut and try solution on a cede plot 



shewed that a gain of 3 X 10 6 and two phase lead 
compensators with a double zero at Z = 70 and a double pcle 
at P = 7CC will satisfy the closed loop magnitude and phase 
requirements. Figure 17 shows the system response and 
desired response obtained for these values. The magnitude 
of the compensated system's response is 93 % of the desired 
response and lags by 8.12°. 
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flGUBE 17 Conventional Lead Compensated System Bespons 



The problem was then run on the optimization program 
with the fclcck connections as shown in Figure 18. The free 
parameters were the poles and zeros of the compensators and 



the gain of the plant. The desired data 


curve 


of 


XBATA 


= sir (200t) was generated 


from the standard s 


econd 


order 


step response by setting 5 = 


0, W = 200 

n 


and 


using 


X (2) 


as the desired response. 


The initial 
/ 


search 


zone 


limits 


were specified as 600 < P 

i 


<800, 60 < 


Z < 

i 


* 

O 

00 


2.6 X 


10 6 < G (3) < 3.2 X 10 6 . 










FIGUFI 



18 CADS Block 



Diagram for Lead Compensation 



The optinizat 



limits 




f cr 


t 


cth 


zeros 


an 


d tt 


e 


gain 


relaxed 


t c 




500 


3.2 X 


10 


* < 


G ( 


2) < 


free 


va 


riat 


le 


s c 



icn program 


solution went to 


the lower 


poles 


and 


the upper limits 


for both 


the 


. The 


limits of the search 


zones 


were 


< 


P 

i 


<600, 80 < Z 

i 


< 9 C 


and 


3.4 X 


10 6 . 


The program again 


placed 


the 


n the 


limits of P = 500, Z 


. = 9C, 


and 



1 l 

G (3) = 3.4 X 1C 6 . The system and desired response for these 
values is shewn in Figure 19. The system magnitude and 
phase are much closer to the desired output response than 
the respcEse of the conventionally designed compensator. 



j 
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FIGURE 19 CADS Lead Compensated System Respcns 



The phase difference is only 4.33°. Continued relaxation of 
the boundaries produced the compensated system shewn in 
Figure 2C. Figure 21 is a plot of the system's response for 
these values. The system's response is improved in that it 
lags the input ty only 3.46° and the magnitude is 
essentially the same as the input magnitude. 




FIGUEE 20 CAES Compensated System Elock Diagram 

The optimization program was activated at T = 0 
although the problem specifications were for the steady 
state response. The problem was rerun with the optimisation 
process started after the transient had died out. The same 
values for the optimum compensator i»ere cTbtained. 
Apparently the snail initial transient did not effect the 
problem solution. 
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FIGURE 21 CADS Relaxed Boundary System Respcns 



E. 



SEE VC S1S1EM COMPENSATION 



The foregoing examples of the simulation - optimization 
program were simple examples which could obviously be 
guickly solved with the standard cut and try methods. A 
mere complex and challenging example cf program usage is 
presented by the system shown in Figure 22. 

The system is an operational servo drive mechanism with 
multiple inputs and discrete level feedback. This highly 
nonlinear system's output was to follow the input as closely 
as possible and in steady state (t> 75ms) there was to be 
very little ccise ripple. The free parameters of the system 
are the poles and zeros of the compensator, C , and the 

Xi 

poles cf the noise suppressor C . 

n 

To simulate and optimize the system it was redrawn using 

the available program blocks as shown in Figure 23. luring 

simulaticr it was found that the current limiter for I was 

D 

net needed because |I j < 25 A and the limiter was removed 

D 

to decrease program run time. The desired response (XDATA) 
was written as a set of three eguatiens. These eguations 
were then used to replace the second order step response 
eguatiens in the standard program. DEVIN (1), (4) and (5) 

were also written as a set of eguatiens and placed in 
subroutine EIANT. The discrete level feedback to blcck two 
was achieved by making DRVIN<2) = INTGEE (THA (12) ) . These 
changes to the main program and PLANT were all that were 
neccessary to simulate this system. The implementation of 
these changes to the program is shewn on pages 80 and 81. 
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JJGUflE 22 Servo Drive Mechanism, Original Compensation 






s 
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FIGURE 23 CADS Dagram of Servo Drive Mechanis 



The system optimization was initially broken into two 
parts. Ibis was to reduce the number of variables that were 
to be optimized per run and to obtain near optimum starting 
values for the free parameters. The above steps were taken 
in an effort to decrease the computer time required for a 
solution. The first run was to optimize the compensator. 



, independent 


of 


the 


noise 


suppressor. The second 


optimization run 


was 


to 


select 


values for the noise 


suppressor, C . 


The 


rationale 


behind this separation was 



that the two circuits perform different functions and should 

therefore be initially separable in their effects on the 
system. The final optimization run was to be made with all 
free parameters available to the program for optimization. 
The search zone centered on the values found above. 

Figure 24 shows the response for the system as 

originally compensated. The optimization run to set the 

values for C resulted in Z = 43.5, P = 21.0, 

L * 1 

Z = 47.5, £ = 5S2.C. Figure 25 shows the simulation of 

2 2 

the system using these values. The initial velocity 

overshoot has teen reduced and the average velocity after 
the transient appears more equally distributed above and 
below the desired velocity. 
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FIGURE 24 Original Servo System Fesponse 



% 
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PIGUfiE 25 CADS Compensated C Systen Eesponse 



The initial optimization run for setting the values of 
the noise suppressor used the same cost function which had 



J ~ J Err 2 dt. This resulted in <f = 0.2 (lower bound), and 
H = 15C0 (upper hound) . Figure 26 shews the result of 



using this cost function was to further decrease the 

overshoot. Ecwever, it allowed larger switching cr noise 
transients as a result of weighting large transient errors 
mere than the lesser noise jitter. Tc overcome this 



This was tc weight the steady state errors mere heavily than 

the transient errors. Using this cost function the values 

set by the optimization program were <f = . 2 and 8 = 1430. 

n 

The damping factor, S , was again placed on the lower limit. 

The initial overshoot was still improved ever the original 
system*s response tut the switching transients were not 
improved ever the previous optimization trial. The results 
of the simulation run using these values is shown in Figure 
27. 



been used for all previous optimization runs 



problem, the cost function 
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FIGORK 26 CADS Compensated C -1 System Respons 



% 
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FIGURE 27 CADS Compensated C -2 System Eespons 



A final optimization run was made with all six 
variables free. The performance index was changed tc the 
weighted cost function shown in eguation V. The weighting 
factor was computed to equally weight the start up transient 
and the steady state responses in an effort to reduce rhe 
switching transients. 



The optimcir compensator values established by CADS were 
P = 21.1, Z = 52.5, P = 521., Z = 48.1, 8 = 0.32 and 

2 2 3 3 

K = 1268. Figure 28 shows the fully compensated system's 
n 

response. The initial overshoot has been reduced and the 

average velocity in steady state is closer tc the desired 
value than the original system's response. However, the 
transient error cn shut off is still present. Table II. 
summarizes the results of the optimization of the servo 
system . 

The CPC time reguired for the optimization process was 
net decreased by splitting the problem into two separate 
parts. The time reguired for each optimization run on the 
individual compensators was the same as the time reguired tc 
optimize the complete system with six variables. 




0.0 < T < 0.045 
0.045 < t < TF 



(V) 
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FIGURE 28 CAES Optimized System Besponse 



Compensator Cost function 
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This problem has presented an excellent example cf the 
necessity tc carefully select the cost function which will 
measure the system's performance. Defining a performance 
index which will weight the more objectionable 
characteristics of a system’s response so that they are 
eliminated or reduced becomes difficult as the system's 
complexity increases. The performance index, J =/ (E rr ) 2 dt, 
was not adeguate in its treatment of the ncise and switching 
transients fcr the above problem. The ether performance 
indexes used were also marginal in their effects upon 
parameter optimization. Although the cost function was 
reduced to a mathematically correct minimum, the system's 
performance was not the best that could be achieved. The 
system's performance was only optimum due tc the definition 
of the cost function. The system was simulated using the 
compensator values determined from the last optimization run 
with the exception of S which was increased to <f = 0.5. 
This simulation was made based upon engineering judgment of 
the effects cf varying £ . Figure 29 shows that this change 
in the damping factor produced a system response which was 
nearer the desired response than any of the previous runs. 
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HGURE 29 Servo System Response, CADS Values with S 



IV. 



DISCUSSION . CONCLUSIONS AND R 5CCMMEND A 110 NS 



A. DISCUSSION 



The objective of this thesis was to investigate the 
feasibility cf developing a computer program which would 
optimize a variety of control systems* responses in the time 
domain. Ihe program developed during the investigation 
proves that time domain optimization cf control system 
responses is Ect only feasible but desirable due tc the 
readily in ter pretable results of the optimization process. 
CADS reguires approximately 200K bytes cf computer core when 
the high speed printer plot is used for graphical display of 
the output arc 230K when outputting calcomp graphs. These 
core reguirements are a maximum and could be reduced by more 
careful, professional programming. 

The computer time required for CADS to arrive at a 
solution is dependent upon 

(1) the order of the system being simulated, 

(2) the area cf the search zone determined by the 
upper and lower bounds on the variable parameters 
and 

(3) the nearness of the starting guess to the 
optimum parameter values. 

Every trial value cf a parameter selected by EOXPLX reguires 
a complete simulation of the system in order to evaluate the 
system's response and compare it with the desired response. 
The program run time is therefore, T Rut4 = the number of 
trials X system simulation time. A high order system may 
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require t*«Dty seconds of CPU time for simulation. If 300 
trials are required to determine the optimum parameter 
values, the total CPU time would be 100 minutes. 



Example 


Lower 


Starting 


Upper 


Number 


CPU Time 


Problem 


Bound 


Guess 


Bound 


of 


Required 


III 


XL 


XS 


XU 


Trials 




b 


-300 


-34.3 


-3 


55 


7min03sec 


(pseudo) 














.001 


.01 


0.1 


225 


24min53sec 


c 


.01 


.1 


1.0 


225 


24min53sec 




400 


450 


500 








400 


450 


500 






d 


90 


95 


100 


2,000 


68min 




90 


95 


100 


• 






3.4xl0 6 


3 . 5xl0 6 


3 . 5x10® 








10 


20 


30 








40 


50 


60 






e 


400 


500 


600 


230 


4hr 




.3 


.5 


.55 








1150 


1300 


1400 










Table III 






Table 


III presents a 


summary of 


the search zones and 


times reguired to obtain a solution for 


seme of 


the problems 


considered 


in this 


thesis . 


The amount 


of CPU 


time reguired 



for a solution, especially for problem III-E, may seem 
excessive. However, there are several considerations which 
should be mace prior tc arriving at this conclusion. 

(1) The equations of the system do not have to be 
written, programmed nor debugged if the system’s 
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component transfer functions are known. 

(2) lime domain requirements do not have to be 
translated irto frequency domain specifications 
for system simulation and design. 

(3) A systematic search is carried cut to obtain 
the optimum parameter settings. This assures that 
with valid bounds on the variables an acceptable 
scluticn will be obtained with the first 
optimization attempt. 

The time reguirea to perform the above steps in the design 
process by conventional means may result ic many more hours 
of CPU time than if the program CADS were used from the 
beginning of the design process. 

Several means of reducing the computer time required 
fcr optimization were previously outlined in section II. 
The most significant reduction is obtained by keeping the 
number of integrations required for system simulation to a 
mirimum. Elcck diagram reduction of the system should be 
accomplished whenever possible. The example of the Sard 
Leonard drive system shewn on page 26 can be reduced to the 
simple system shown in Figure 30 by block diagram reduction. 




r 




El CURE 30 Elock Reduced Ward-Lecnard Drive 
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Section III.E. showed the effect of using different 
cost functions to measure the quality cf the optimized 
response. Although the integral error squared criteria is 
often used tc judge a system's performance, the user should 
carefully ccrsider how to best define a cost function which 
will properly penalize deviations of the system response 
from the desired response. A small error (err < 1) squared 
becomes even smaller. If all the errors are small, the IES 
is a valid cost function but if the system response involves 
large and small errors a weighted cost function will have to 
be used. One method cf arriving at a properly weighted cost 
function is tc simulate the system using first estimates of 
the variable parameters and recording the sum of the errors 
over the different portions of the response. A ratio cf the 
errors can then be used tc arrive at proper weighting 
factors fcr each time section cf the response. 

BOXPIX will continue the optimization process, working 
in the seventh significant digit until it can no longer 
reduce the cost function. Often an acceptable solution for 
the system parameters has been found lorg before the cost 
function has teen reduced to its minimum. A judicious use 
cf CADS nay be made by evaluating the system response after 
ten to fifteen minutes of run time to see if an acceptable 
scluticn has been found. 



B. CONCH SICKS 

lime doaain cptinization using the CADS program is a 
simple, straight forward process which dees not reguire an 
in-depth analysis cf the system being optimized. Economic 
use of CPU times does dictate that intelligent starting 
values and bounds be placed on the variable parameters which 
are to be optinized. 
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The program is a readily usable tool for simulation 
without cptinization. It is competitive with, if not 
superior to, other common simulation routines when 
simulating typical control systems. This feature alcne is 
expected to bring the program into common usage. 

C. RECCfiMENI ATIO KS FOR FUTURE WORK 

(1) All integral calculations performed during this 
investigation were done in double precision for increased 
accuracy cf the system response. The possibility of using 
single precision calculations should be investigated to 
decrease core reguirements. 

(2) The ability to begin the optimization process at some 
time greater than zero should be provided. This will 
necessitate modification of the data input cards and block 
eguations sc that initial conditions can be entered. 

(3) The program presently reguires that external forcing 
functions (ERVINs) and the variables which are to be 
optimized be specified by Fortran IV statements placed 
within the bcdy cf the program. A method cf reading these 
specificaticrs frcm data card input should be developed sc 
that the user will not have to ’’shuffle" cards in the 
pregram deck. 

(4) The graphical output of CADS was all that was necessary 
rcr the investigations conducted in this thesis tut a 
prevision fci numerical output of selected responses should 
be -provided for detailed analysis. 

(5) The feasibility of reducing the number of significant 
digits ECIFIX uses should be studied as a means of reducing 
optimization time. Also some criteria might be developed to 
reject a system' s response before TF is reached if it is 
determined tc he unacceptable. 
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(6) A method of automatically relaxing the boundaries on the 
variables being optimized when they go to their limits 
should be developed. 

(7) The standard cost function provided and all user 
developed cost functions should be normalized. This would 
permit a more direct and easier comparison of a system's 
"goodness" when several different integration step sizes or 
run times have been used in the optimization process. 
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Nirricn 



APPENDIX A 



Block Data Card Format 



111 1 20 1 40 1 6 0 



ELKCC D= EE V V 



G 



P 



Z 



CC = POSITION NUMBER OF THE BLOCK 
E = TYPE OF ELOCK (NUMBER) 

EE = INPCI NODE NUMBER 
V V = OU1EDT NODE NUMBER 

= VAIOE OF GAIN 

= VAIOE CF THE POLE* 

= VAICE CF THE ZERO* 



* Fcr fclcck type 4, S is read into the P locaticn and 
Wn is read into the Z position. 



* Ecr tlock type 5, 0 is read into the P positior 

U 

and 0 is read into the Z position. 



* Ecr ticck type 6, © • is read into the P position 



and Q e is read into the Z position. 



EXAMPLES 



G 



El KC 11 = 0 102 



10 



G 



P 



BIKC22=03C4 



1 . 5 



G 



P 



Z 



BI F C23=0 405 



10. 5 
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El K 1C4 = 1 C06 



BIK1 15=1207 



El K 1 26=0 7 08 



G b u), 

200. .2 14.14 

G 0 U ~© L 

10. 20. -5. 

G e c e fi 

1. 3. 10. 



♦See Table I. 
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APPENDIX B 



Plant Flow Chart 



FUNCTION PLANT (C) 



CINENSICN G<2 5), P<25), Z(25>, FLAG(25,25), CM 
C-(25), TFACC7 ( 25 ) , G N GCOT (251, DPVINM25), 

N F ( 2 5 ) » IP (25) » IV ( 25 ) ♦ X2(25), X2DOT(25)t 
CFIVE(25), TFAI25), IC(25>, IC(25), IE(25), 

TFACLT (2CC1 ) , XC AT A ( 300 1 ) » C(25) 

REAL *87 F ACL7 

REAL * € X C A T A 

REAL *87FA,TFAC07,T ,DT,CRIVE ,C N G , C NC-CGT , CR V I N ,TF 
REAL *8F1,F2 
REAL *SX2,X2CCT 

CCNMCN 7 ,CT,7F ,7FACUT , XDA7 A , M 3 , I CCNT ,NECt ISK IF , IIP 
FCP OPT RUNS, IN FI E IT REAC STATEMENTS AFTER REACING 



* 

* 



* 

* 



* * • 

IF * 
ISK IP-1 



* ? * 
0 I 



I 1 



I 



. * 

+ 

T 5 



1 I ISK IP = 2 



I 

*** F E A C ( 5 , 2 A ) N » IS ET 
INITIALIZE CCLNTEFS 



I CUT =C 
I \CLT=C 
F 1 = CT 

F 2 =C.5CC*F1 

Ml = C 
N 55 = C 

N66 =0 

ICK = 2 * N 



I 
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FEAC IN FIT CATA 



4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 



4 — « 

4 CC 

•2 

4 i«i;n 

4 



I 

***REAC (5,25) IC(I)fID(I)fIE(I)flV(I)tG(I)fP(I)f 
SET TFACIT = C IT FIT CF LAST ELGCK 



Z(I) 



• * * • 

• * IF * . 

*(l\( I J-IVCLT ) • LE«0 * 

^ ■ 

* . . * 

* . * 

wm 



nciT=iv(i ) 

ICIT = 1C <! ) 



* 

* 



• * ' * • 

* IF * 

ica > . ec • i 



* 

*- 



4 

* 



IF 



* 

* 



$ 

* 



IF 



I 2 



I M1=N1 1+1 



* 

1 


1 


* . 

* 3 

EC. 5 * 


. * T | N 5 5 = 


=N55+1 I 


1 


1 


* . 

* . 

cC.6 * 

» ^ 

. * T | N66 = 

4 fc — 


= N£6 + 1 | 



4+4 + 4 , 44^R1TE (6,26) ICm,IC(I),IEm,IV(I),GC>,FlI),Zm 
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SET TFACLT 



= SPECIFIED ThETA, IF ANY 



IF 

ISET.NE.C 



*- 



***WRITE ( 6 t ^ 7 ) I CUT 



Ml = N 1 1 
^ 5 5 =^^N5 C 

N66 =4*N6C 

NEC - N - 1 



I 



SCAN FC P INFLTS 



I IGLT = 1 SET 



CONNECT SYSTEM E\ SETTING FLAG = 1 • FC F CCNNECTEC ELKS 

I 

4 + 

« cc 4 

+ + ^ 4 4 

4 + J - 1 t N 4 

4 4 4 

4 
4 

4 4 

4 4 CC 

44444+4444 4 



4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 4 
4 

4 

4 

4 

4 

4 

4 

4 



4 K = 1 t N 
4 



I F L AG ( K f J ) =C.O 



* IF * 

IV(K).EC.IE(J) 



* * 

: 



* 

T 



I FLAG ( K 9 J ) = 1 • C 



TAKE eLK TYPE ANC SCLVE ECNS ONE EY CNE 
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CLEARING CLT F EG ISTERS AN C INITIALIZING COUNTERS 



4 

♦ 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

t 44444444 



4 

4 ‘ t£ 

6 

4 ICLP=1»N 

4 



ThA(ICLR) =C.CO 
ThACCT ( ICLF ) =0 • CO 
CMG(ICLP) = C.CO 
CMC-CCT (ICLF ) =c.co 
CFVIN(ICLF) =0.C0 
X 2 ( I CLP ) =C.CO 



>2CCT ( 1CLP ) =O.CG 
T =C.OCC 



NR'S CCNTFCL ENTRY FCIM IN INTEC-RATICN SUEPOUTINES 



NF2 =1 
N P 2 = 1 

N F A = 1 



N'2 CCN7FGLS V- 1- 1 C h BLOCK EOLATION IS BEING SCLVEC 



I M 2 =0 



ICCNT IS LSEC TC CCNTFCL FRCGRAM FLOW 



I I CCNT = C 



IWA IT IS LSEC TC CCNTRCL TIME. TIME IS STEPPEC EVERY 



FIRST £ N C TMPC FASS THRU INTEGRATION ROUTINES 



I WA IT = C 
IT =C 



I LAST 1 S CCNTFCL FFCGRAM CATA AND PROGRAM EXIT 



I L AS T = C 
1 SL AST =C 
I 6 LA ST =0 
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SPECIFY AN > F^FAFETEFS TC BE OPTIMIZEC HERE 

i 

EG • F ( 1 J=C(2), Z(1)=C<1), 

SET CRIMES FC? 1NFITS 



I 

4 4 

* cc 4 

-♦4444 + 4+44 <5 4 

4 4 PCRVMtN + 

4 4 4 

4 
4 

4 

4 I CRVIN(l) =1. 

4 I CFIVE(KFV ) =C.CC 

4 

4 
4 

4 4 4 

4 4 CC 4 

444 i *444444 8 4 

4 4 f* = 1 t N 4 

4 4 4 

i 
4 

4 

44 4 + 4444) CFIVE(NCRV) =CRI\5(NCRV)4TM(M) 

4 j *FLAG(N*VCFV) 

4 

•» I 

•» 

44444444} CP I \i 5 ( NCR V ) =CRVIN(NDRV)+CRIVE(MCRV) 



10 I N2 = N2+1 

- 

PICK 7 \ F E ECN TC SCLVE 



* 

5} 



* 

* 



* IF 
IfcMT.EC.C 



* 



• * 



4 IF * 

IC <N2 J.EC.l 

* . • * 

* ! 4 
F 

! 



if 

if- 



T=T+H2 



I 11 I 
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-tf- -W- -W- 4i. 41. -H- -W- 



• * 4 . 

• * IF * . 

I C ( V 2 ) . E C . 2 * 

• • ^ 

* • . * T I 12 

* 1 * 



. * IF * . 

IC(N2).EC.2 * 

• • ———————— 

* . . * T I 12 

* 1 * 

F 



• * * • 

. * IF * . 

IC (N2 1.EC.4 * 

• • $ — — — — — — — — — — — — — — 

* • . * T | 14 

=> . * 

I 

F 



. * IF * . 

IC(N2).EC«5 * 

a • >!< — — — — — — ” — 

* . . * T I 15 

^ # ^ — — — — 
I 

F 



• * * • 

. * IF * . 

IC(M2).EC.e * 

• • ^— — — — — ——————— 

* . . * T I 16 

^ Q 5 



***WRITE (6,28) 



S7CF 



I 



74 



ST 4 F 7 SCllTICN 



TYPE C N E EC14TICNS 



11 



SOLVES Tl-ACIT = G*THA IN 

I 



7I-A(M2 ) =G (M2 )*CFI VE(M3) 

IWAITMMIT-H 

ILAST=ILAS7-U 



• * ' * • 

IF 4 

*(flAST.EC.Nll).ANC.(ICCM.5C.NEQ) 

* • • * 

* . . * T |2 

*-» . * 



I 18 I 



12 



12 



TYPE TkC ECNS 



SOLVES ThAClT/TPAIN = C-/(S + P) 



TMACGT (M2) =-P(M2)*ThA(M2)+G<M2) 
*CPIVS (M3 ) 

S =FKLCE2(TPA ,TI-A00T,NR2 ) 
IMIT = IMI7 + 1 



• * * • 

. * IF * . 

* S-l.C * 

* . • * 

* . 4 



I 17 I I IE I I 21 j 



TYPE 7PFEE ECNS 



SOLVES ThACLT/TMIN = G*(S + Z)/(S + F) 



CMGCGT (M2) »-P (M2 ) *CMG ( M3 ) +C- (M2 ) 
4CFIVE (M3) 
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I 

$ = FKLDE2 (CNC- ,CNC-DCT,NR2 ) 



TFA(M2> =<Z (K3 )-F(N2) )*0MG<N2) 
■*G(N3 )*CPI VS(f*3 ) 
UAIT=IVAI7-U 



* 

* . 



I_17__l 

1 



. * * 
IF 

S-l.C 




I 16 I 



* 

* 



* 



I 21 | 



7V F E FCIP EONS • 

SOLVES THACLT/TFAIN = C-/(S**2 + 2*DEL7A*WMS + kN**2> 



14 



V =CCFLXF ,Z, C-, CRIVE, X2 ? CNG,NR4 ) 
TFA(F*2) = CK-<N3) 

IV»AIT = lVAn + l 



• * * . 

. * IF * . 

* V- 1 • * 

* . . * 

* • • * 

* . * 

0 I ♦ 



I 17 I 1161 1211 



I 

TYPE FIVE ECJvS 



SOLVES TFAQLT = G*ThMN FOR /G^TFA IN/ < SAT LIMTS 



15 



| T h A ( N2 ) »G(P3)*DPIVE(M3) 



I 



• * * • 

. * IF * • 

* TFA(N2)«CT.F(f l '2) * 

^ * * . . * T T TH A ( M3 ) = P ( N 2 ) 

4 . 4 

I 

F 
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a =» * • 

. * IF * . 

* T F A (M2 I • LT . 2 (M3 ) * 

^ a a $ — — — — — 

* • . * T | T F A ( M 3 ) = Z ( M 2 ) | 

* , * 

I 

F 



I I l*A IT = IV»A IT + 1 
I I £1 AST =I5LAST + 1 



. * * . 

• * IF * . 

* ( I5LAST.EC.N55 ). ANCa ( ICGM.EC.NEG ) 

^ • • 5 ^ — — — — ——— — — — — ————— — — 

* . . * T | 21 | 

* , * 

F 



I 1£ I 



TYPE SIX ECNS 



SCLVES T FAG IT = G*TFA I N FCR CGNT REF > C Z 3CLNC 

I 



16 I TF A ( M 2 ) = G (M2 )*CFIV£(M3) 

1 IF = F ( N 2 ) 



. * * . 

• * IF * a 

*( C A ES ( TFA ( IF ) ) ).LT.Z(M3) 

i a 3^ — — — — — — — — — — — — — — — — — — — — — — 

=» a a * T | TFA ( M3 )=0a 

* . * 

/ I 



I I lnAIT = IfoA IT+ 1 
I I6LAST =I61AST+1 
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• 4 4 • 

. 4 IF * . 

* (I61AST.EC.N66 ).ANC. ( I COM. SC .NEC ) 
4 . .4 — — 

* . . * T I 21 

4.4 

I 

F 



I 16 I 



17 



16 



19 



*44^PITE (6,29) 
STCF 



I ICCNT= ICCM-fl 



4 

4 



. * * . 

* IF 4 

ICCNT-r 



4 

0 



4 

4 



I 1C I 



19 

T 



20 



CNE PASS HAS E £ E N PACE FCP EACH ECN. 

NP2 = N P 2 + 1 
I' P 2 =NF2+1 

NF4 = N P -A + 1 
P 2 = C 

I CCNT = 0 



INCREMENT 



4 

4 



• 4 4 . 

4 IF 4 

IVnAIT.EC. I C K 



* 1 * 
F 



• 4 



1 7 I 



I I WAI T=C 



2C 



444ViRI TE (6,20) 
ST C F 
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21 

22 

25 

24 

25 

26 

27 

28 

29 

2G 

EN 



CATA CCLLECT ICN FCINT 



CAT A IS RECORDED AT 5NC IF 



COMPLETE TINE STEF(CT) 



I IT =17+1 
I TFACUT ( IT ) =TH A ( ICLT ) 



I 

TEST FCR ENC CF FIN 



. * IF * . 

* T-TF * 

^ . • * 

* • • ^ 

* . * 



I 22 I I 22 I I 25 | 



I 

RESET CCLNTEFS FCR NEXT PASS THRU EOLATIONS 

I 



N R 2 =1 

NR5 =1 
NR4 =1 
I* 5 =C 
I CCNT = G 
I W A IT = C 
I L A ST = C 
I 5L AST =0 
I6LAST =0 



I F L ANT= 1 • 



RETURN 

FCFNAT ( 1 2 » 2 > 1 12 ) 

FCFNAT ( 5 X » I 2 ♦ 1 1 » 1 X » 2 1 2 » 8X » 2 E 20 ® 7 ) _ 

FORMAT (/»5F ELK »I2fIl»2H = » 21 2 » 6X» 3E2C. 7 ) 

FCFMAT (/ / »2X» ' T F £ 7 A CUT IS TH A { * , 1 2 , • ) • ) 

FO F NAT (4CH *4* EON SWITCH CONTROL DID NCT WC^K) 

FCFMAT (2CF INTEGRATION TRCUELE) 

F r FMAT ( C EF ERROR IN INTER CATION. ATTEMPTED TO 
IN7EG. MCFE THAN N-EGN) 
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APPENDIX C 



oooooooqoooo 
O — < fvj m ^ vo oo cn o «-* 
Cj>U>gr>^nu>cr*cr*Lr'U > U''00 



OOOO 
CTO — «<\J 

cvimmoi 



-a 

t— 

c 

3 

uo 



LL 

X 

o 

LU 



cr 

a 



iu 

cn 

2 

3 
LL 
00 



CL) 

o 

a_ 

CL 



UJ 

CL 







2 . 


















2 








UJ> 




*— 4 














4-H 








H- CL 












UJ * 






H 








l/)3 




II 








3 














>- 












2 






U. 








l/) + 




LL 








M 






3 




0 








C 








H 










UJ 




QC'X 




UJ 








2 






LO 




1 




LU-fr 




3 








LJ 






2 




3 




C-J CM 












LJ 






O 




4—1 




ai-fr 




LL 














4—4 




M 




O* 




1-4 














h- 








2 








H 










o 


<r 


X 




U3 




LU 




r -4 










2 


c\J 


UJ 




2 




;> 




0 ) 00 




U) 






3 


l\i 3 


3 




3 1 




Cti 




HCVJG 


— * 


UJ 






LL 


c\Jn)H 


aj 




LJ 




3 


in 


• 00 3 


H 


3 


H- 






•—4 <\J 


O 


2 


LUO 




o 




— <C0O 


• 


3 


2 




3 




QC 


3 


oox 






u 


« M" 


O 


< 


C 




4—4 


J>X -4|— 


CL 


C 


* 




c 


»— 


o 


| 


m j> 


3 






cr +3 




X 


02 




1— 




i—ao 


1— 


ao 


CL 




3 


UHHUOI- 


CC 




iT2 




< 


o 


H-O 


*— 


—••—4 0 






CC 


2 II H-O** 


o 


a 


<*fr 




o 


CD 


CJ 


*■ 


»— l »*oc 


O 




IU 


-'*-•1-0 O 


LL 


h» 


Gu 3 








oa— 


0 


— m< 


1- 




X 


OJ— OO — vf 






2* 




UJ 


— * 


OOh 


O 


xcua 








— iH-— * OO • 


0 


oo 


CCNJ 




1“ 


r-t 


— m* 


s0 


-4 II -*2 


V) 


2 


o 


— <*dO— ' rOH- 


UJ 


2 


h- 1 




3 


• 


cn— •oj'U 


o 


+ — **C 


2 


► 


CD 


CN 4 XO •O — C\KT 


3 


a 


</) 


O 


0 l 


30 


330 ^ 


t-4 


CC3|— 


LJ 


3 




Hf- JU^GQOI 


LL 




1 ! 


o 


X 


LUO 


• • • •'T 


<T 1 >T 


H-^-cuoo 


M 


II 


cn 


33 O.UJ • • ••— 4 


M 


H- 


LU 


o • 


O 


• • 


OOO^^cohoj^ 


<C3 


I— 


> 


— 


Cl— H »ooo H 


LJ 


C 


f— LJ 


o 


LJ 


CLLJ 


• • •'TOO 


• uj «aj 


CJCJ — 'UJUJ 


C 


UC 


— 


-L2-*H- • • ' 


UJ 


O 


<LJ 


II 




< 


UJUJ LL) m «— 1 




• HWU.>0 


LJ 


3 


4—4 




2. 


*•«* 


GX 


11 


H- 


UJ || OH- 030-4 


m m 


O H — H<tZ 


*— < 


2 


— 


II II — C333— 




LL 


3 


c 


O 


3 


CJ t • • II O ll a II O II c<i 1 X«~* 


LL 




2 


M>23 e • »2 


2 L 


M 


LJ 


I-—* 


2 




ll thhh^h 


— H— *H- 


— Hh*H H" 


•—4 


U> 


3 


•— (X »-«* 


UJ 


LJ 


3 


< 1 — « 




IM 




*-H — 1 




3 




-> 




UJ 


LJ 


C 


3^ 


LJ 


LL'— 1 


II U.U-U-— u 


LJ— LJ 


'-'uuu-uU 


3 


3 


Lt_ 


X^U-U-U-LLLta. 


LJ 


2. 


LJ 


3X 


CJ 


wX^h’-’wwXUXUX^XMXw^U 


2. 


3 


3 





au 

r-^ 

uuuuuuuu uuu 



r— 4 


IM 


t*»sr 


01 






3 


3« >4 


UJ 


3 


au 


3 


00 00 


3 


3 


— i 




33 




3 



o uuuu ooo 



80 



oooooooo 

ro in >o r- oo o 
mcnrorororom>r 



oo 



f-l 


vO 






O 




o 






• 


CM 






ro 




W 






o 






oo 


«-H 




ro -ft- 






1 






CM 


hh 




0)00 






*— 


CM 






ro 




nr- 






w 






H 


o H 




— KJ> 






-fr 


1 — — * 






hH 




r- 






• 






o 






O • 






o 


uo 




1 — 


OLJ 


f— 


1 — CM 




o 


o 


Oh 






C 2 h- 




ro 




a 


o 






a 




CM 


o + 




• 


rH 


— o 




OCT* 


— 00 s 


CM 


OO 




o 


1 


vOO 


O' 


o 


vOOrO 


cn 


ro 






00 


ro 


o 


—co 


ro oo 


v 0 


—CO 




ii 


• o 


CM~» 


00 


ro • 


(M — CT 


• 


rocr* 






H • 


OH 


• 


f-HCM 


• OH « 


i-H 


r—4 • 


• 





tnrotnromo 

c\j u cm n CM !l 



• •OCM'TCM «0(MO 
oo ii f o n roo n o> ii 

J • • — r-M— r-M • — «— I— U 



> 

o' 



• • COCM(J> • I CTO 

OOCNJ II CMO H CM II 

j § »i .< — H 

«in< 

W W W 2 JO w O'— <— O'— 2 T — 

020Z02H • »OZO ®ZU2H • ^ZaZO *ZOZH UJ 

-> .> p-z— > p- — m 

UU.UUiUU.UU.u.uCLUU.U.UU.U^-^^UCX.UU.U.UU.U UL 
UUOUOUUHHOUUHUOU(J*-'HUOUU»-»UOUU lj 



II 



<JU 



nT 

CM 



IM 

iM 



1 NU 1 

CVJCM 



vO 

CM 



UJLP 

tMlVJ 

HH 



o 

PI 



IMIV 

1 MO) 



o o 



o o 



81 



CRIVE(MCRV) = DRIVE(MORV)+THA(y ) * F L AG ( M , MC P V ) 
CFIVE(MCRV) = ORVIN(MCRV)+DRIVE(RCRVJ 
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KE EVALUATES IMPLICIT CONSTRAINTS ON BQXPLX VARIABLES 
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